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Abstract 

Genetic diversity generally underpins population resilience and persistence. Re- 
ductions in population size and absence of gene flow can lead to reductions in 
genetic diversity, reproductive fitness, and a limited ability to adapt to environmen- 
tal change increasing the risk of extinction. Island populations are typically small 
and isolated, and as a result, inbreeding and reduced genetic diversity elevate their 
extinction risk. Two island populations of the platypus, Ornithorhynchus anatinus, 
exist; a naturally occurring population on King Island in Bass Strait and a recently 
introduced population on Kangaroo Island off the coast of South Australia. Here 
we assessed the genetic diversity within these two island populations and contrasted 
these patterns with genetic diversity estimates in areas from which the populations 
are likely to have been founded. On Kangaroo Island, we also modeled live capture 
data to determine estimates of population size. Levels of genetic diversity in King 
Island platypuses are perilously low, with eight of 13 microsatellite loci fixed, likely 
reflecting their small population size and prolonged isolation. Estimates of heterozy- 
gosity detected by microsatellites (H E = 0.032) are among the lowest level of genetic 
diversity recorded by this method in a naturally outbreeding vertebrate population. 
In contrast, estimates of genetic diversity on Kangaroo Island are somewhat higher. 
However, estimates of small population size and the limited founders combined 
with genetic isolation are likely to lead to further losses of genetic diversity through 
time for the Kangaroo Island platypus population. Implications for the future of 
these and similarly isolated or genetically depauperate populations are discussed. 



Introduction 

Genetic diversity has been identified as an important fac- 
tor influencing a population's long-term potential for sur- 
vival (Bouzat 2010). The contribution of genetic diversity has 
been recognized in numerous aspects of population persis- 
tence, and is critical for long-term fitness and adaptation [see 
Frankham (2005) for a review] . A loss of genetic diversity has 
been shown to affect individual fitness with decreased sperm 
quality (Hedrick and Fredrickson 2010), reduced litter size 
(Hedrick and Fredrickson 2010), increased juvenile mortal- 
ity (Ralls et al. 1988), and increased susceptibility to disease 



and parasites (Coltman et al. 1999). Accordingly, popula- 
tions lacking genetic diversity often exhibit an increased rate 
of extinction (Markert et al. 2010). Inbreeding, genetic drift, 
restricted gene flow, and small population size all contribute 
to a reduction in genetic diversity. Fragmented and threat- 
ened populations are typically exposed to these conditions, 
which is likely to increase their risk of extinction (Saccheri 
et al. 1998; Madsen et al. 1999; Frankham et al. 2010). 

Island populations are often isolated and small in size, 
and therefore experience increased levels of inbreeding and 
a greater impact of genetic drift. As a result, these popu- 
lations generally have lower levels of genetic diversity and 
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fitness than counterparts found in mainland populations. 
For example, the black-footed rock- wallabies of Barrow Is- 
land, Western Australia, have extremely low levels of genetic 
diversity (H E = 0.053), which has led to this population 
suffering inbreeding depression: females experience reduced 
fecundity and individuals exhibit increased levels of fluctuat- 
ing asymmetry (Eldridge et al. 1999). These findings are not 
atypical with island populations frequently identified as hav- 
ing reduced genetic diversity, increased levels of inbreeding, 
and therefore a higher extinction risk (Frankham 1997). 

The platypus, Ornithorhynchus anatinus, is found in east- 
ern Australia, from as far north as Cooktown in Queensland, 
along the east coast to Victoria and throughout Tasmania 
(Fig. 1). Two island populations of O. anatinus exist within 
this distribution; King Island in Bass Strait, and Kangaroo Is- 
land off the coast of South Australia. The population on King 



Island is naturally occurring and likely dates back to the exis- 
tence of a land bridge across Bass Strait connecting Tasmania 
with the mainland. Intermittent rises and falls in sea levels 
have seen this region periodically disconnect and reconnect 
with the mainland until its final isolation about 1 1,800 years 
ago (Blom 1988). King Island remained connected to Tas- 
mania until approximately 10,000 years ago (Blom 1988). 
Despite its more recent connection with Tasmania, prelim- 
inary microsatellite analysis indicate platypuses to be more 
closely related to Victorian individuals (Akiyama 1998). 

The Kangaroo Island platypus population was established 
with introductions from both Tasmanian and Victorian in- 
dividuals. In 1928, one female and two males were trans- 
ferred from Wynyard, Tasmania (Anonymous 1934; Flora 
and Fauna Board of South Australia Report, 1928-1929). 
This was followed by several introductions from Healesville 
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in Victoria with five males and five females relocated in 1941 
(Anonymous 1941; Fleay 1941; Anonymous 1946) and three 
males and three females in 1946 (Anonymous 1946). Today, 
Kangaroo Island is home to a self-sustaining population of 
platypuses. Akiyama (1998) conducted a preliminary genetic 
investigation of Kangaroo Island platypuses using microsatel- 
lite markers. Genetic diversity appeared to be relatively high, 
but the absence of comparable analyses for other platypus 
populations makes the interpretation of these results difficult. 
The contribution of founding individuals from both Tasma- 
nia and Victoria remains contentious. Previous studies have 
revealed mitochondrial DNA (mtDNA) and microsatellite 
markers to be more similar to Victorian individuals (Gem- 
mell 1994; Akiyama 1998) but, critically, comparisons were 
made to individuals outside the regions from which founders 
were sourced. Current levels of genetic diversity in the island 
population and the persistence of Tasmanian and Victorian 
genes have yet to be adequately addressed. 

Kangaroo Island imposes a restricted area of suitable 
aquatic habitat and isolation from other platypus popula- 
tions. As a consequence, small population size and a lack of 
naturally occurring gene flow are likely to lead to inbreed- 
ing, the fixation of alleles, and associated reductions in ge- 
netic diversity over time. In addition, the limited founder 
number is likely to contribute to problems associated with 
inbreeding. Distribution records on Kangaroo Island suggest 
that platypuses are confined to the western end of the island 
(Fig. 1), further restricting population growth. Knowledge of 
the population size on the island can further inform estimates 
of genetic diversity and give some indication of population 
fitness. Population size relates to genetic variation (Frankham 
1996) and is also significantly correlated with fitness (Reed 
and Frankham 2003). Consequently, small population size 
can reduce the evolutionary potential of a species (Frankham 
1996). Additionally, effective population size can be inferred 
from population size estimates (Caballero 1994) and be used 
to predict the loss of genetic variation over time (Frankham 
1995b). Due to the limitation of suitable habitat on Kanga- 
roo Island, population size is likely to be small and associated 
threats to genetic variation and population fitness are antici- 
pated. 

Here we aim to investigate current levels of genetic diversity 
on King Island and Kangaroo Island by genotyping individ- 
uals at 13 polymorphic microsatellite markers described in 
Furlan et al. (2010). We also sequence fragments of the mito- 
chondrial genes cytochrome oxidase subunit II and cytochrome 
b to determine the relationship between King Island and Tas- 
mania or Victoria and the contributions of founding indi- 
viduals on Kangaroo Island. In addition, live capture data are 
modeled to determine population size estimates for Kanga- 
roo Island and provide insight into the likely future losses of 
genetic diversity in this population. 
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Methods 

Genetic analysis 

Platypus live-capture surveys were conducted using either 
gill nets [refer to Grant and Carrick (1974) for method- 
ologies] or fyke nets. Fyke nets were set in pairs, back to 
back with the cod end secured well above water level. Nets 
were set approximately 2 h before sunset and checked regu- 
larly throughout the night. Individuals were identified by a 
uniquely coded passive integrated transponder tag (Trovan, 
Trovan Ltd., East Yorkshire, UK, or Allflex, Allflex Australia, 
Brisbane, Australia). Trapped individuals were scanned for a 
tag and any untagged individuals had one inserted subcuta- 
neously between the scapulae following protocols in Grant 
and Whittington (1991). Each unique individual had a sec- 
tion of skin webbing approximately 2 mm 2 cut from the 
distal margin of webbing on the rear foot using sharp, ster- 
ile scissors. All samples were immediately placed in 2 mL 
Eppendorf (Hamburg, Germany) tubes containing 100% 
ethanol. 

Platypus web-tissue samples were collected from 12 in- 
dividuals on Kangaroo Island between December 2008 and 
December 2009. Eleven platypuses were sourced from Rocky 
River and one from Breakneck River (Fig. 1), representing 
the only two rivers in which platypus presence has been con- 
firmed on the Island. These reside within the Flinders Chase 
National Park and the adjoining Ravine des Casoars Wilder- 
ness Protection Area, supporting a total land area of 739 km 2 
(Smith 1995). Additional samples were collected from Kan- 
garoo Island (n = 3) and Warrawong Sanctuary (n = 10) 
between 1991 and 1996 as part of a previous study (Akiyama 
1998). Warrawong Sanctuary is located on mainland South 
Australia and represents a captive population sourced from 
Kangaroo Island individuals. Three of the sampled "Warra- 
wong Sanctuary" individuals were originally sourced from 
Kangaroo Island, while the remaining seven individuals were 
subsequent generations bred within the sanctuary. Eighteen 
individuals were sampled from King Island in January 2009 
(Fig. 1). Additional deceased samples (n = 3) were obtained 
opportunistically between 2000 and 2007. The island en- 
compasses around 1098 km 2 of land, although platypuses are 
primarily restricted to the rivers flowing to the east and south 
(Barnes et al. 2003). On the mainland, 60 individuals were 
sourced from the upper Yarra catchment in Victoria (here- 
after referred to as Victoria) between 2008 and 2010 (Fig. 1). 
In Tasmania, 19 individuals were sampled from Duck, Black- 
Detention, Inglis, Cam, and Leven catchments in the states 
northwest (hereafter referred to as Tasmania) between 2007 
and 2009 (Fig. 1). Tasmanian and Victorian individuals have 
been sampled from areas that potentially provided founding 
individuals for Kangaroo Island and have been included here 
for comparative purposes. 
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DNA was extracted from web-tissue using a cetyl trimethy- 
lammonium bromide (CTAB)-phenol/chloroform extrac- 
tion (Sambrook and Russell 2001) and genotyped at 
13 microsatellite loci (Furlan et al. 2010) and two partial 
mitochondrial genes. Reactions contained 2 /xl of template 
DNA, 1 x polymerase reaction buffer, 2 mM MgCl2, 0.2 mM 
deoxyribonucleotide triphosphates (dNTPs), 0.5 mg/mL 
bovine serum albumin (BSA) (New England Biolabs, Ipswich, 
MA), 0.03 U of Taq DNA polymerase (New England Biolabs), 
0.3 /xM forward primer end-labeled with [c33P] adenosine 
triphosphate (ATP), 0.2 /xM unlabeled forward primer, and 
0.5 /xM unlabeled reverse primer. Reactions were made up 
to a final volume of 10 /xl with ddH 2 0 Polymerase chain 
reactions (PCRs) were carried out with an Eppendorf Mas- 
tercycler (Eppendorf, Hamburg, Germany). Samples under- 
went initial denaturation at 95°C for 10 min, followed by 
35 cycles of 95°C for 20 sec, 55°C for 30 sec, and 72° C for 
30 sec. A final extension step of 72°C for 3 min completed 
the reaction. PCR products were then run through 5% poly- 
acrylamide denaturing gels at 65 W for 2-4 h and exposed 
to autoradiography film (OGX, CEA, Sweden) for one to five 
days. Alleles were sized by comparison with Agtl 1 ladders 
(Promega fmoL DNA Cycle Sequencing System, Madison, 
WI). 

Mitochondrial DNA primers were derived from the 
published O. anatinus mitochondrial genome (Warren 
et al. 2008). A subset of samples was analyzed from 
each location. Sequences were amplified to obtain a 
603 bp segment of the cytochrome oxidase subunit II gene 
(COII) and an 835 bp fragment of the cytochrome b gene 
(cytb). PCRs were performed using primer pair com- 
binations COIIF [5'-AATGGCCTAYCCYCTYCAAC-3'], 
COIIR [5'-CCGCAAATTTCTGARCACTG-3'], and 
cytbF [5'-CCCACCCCCTCTAACATCTC-3'], cytbR 
[5'-TAAGGATTGARGCKACAAGG-3'] . 

PCR amplification was performed in 30 /xl reactions con- 
taining 1 x polymerase reaction buffer, 2 mM MgCh;, 0.2 mM 
dNTPs, 0.03 U Taq DNA polymerase (New England Biolabs), 
0.5 /xM of each primer, and 5 /xl of template DNA. Reac- 
tions were made up to a final volume of 30 /xl with ddH 2 0. 
PCR amplification of double-stranded product was carried 
out with an Eppendorf Mastercycler (Eppendorf, Hamburg, 
Germany). PCR conditions comprised of an initial 7 min de- 
naturing step at 94° C, then 40 cycles of denaturation at 94° C 
(20 sec), annealing at 54°C (30 sec), and extension at 72°C 
(60 sec), with a final extension step of 72°C for 1 min pre- 
ceding an indefinite hold period at 4°C. PCR products were 
sent to Macrogen laboratories (Seoul, Korea) for purification 
and sequencing on an ABI3730 XL DNA sequencer (Applied 
Biosystems, Foster City, CA). To ensure accuracy, all ampli- 
fied products were sequenced in both the forward and reverse 
direction using primers from the initial PCR reaction. Of the 
amplified segments, 518 bp of COII and 738 bp of cytb were 



used in the analysis. Negative controls were used to ensure 
contamination was avoided. 

Individuals sampled from Kangaroo Island in December 
2008 and December 2009 were used to assess contemporary 
levels of genetic diversity. The contribution of the founding 
individuals was determined by incorporating all Kangaroo Is- 
land and Warrawong individuals (sampled between 1991 and 
2009). Due to the long-term isolation and low genetic diver- 
sity present on King Island, analyses of this population were 
conducted incorporating all samples (i.e., sampled between 
2000 and 2009). 

Given the genetic similarity among individuals of the upper 
Yarra catchment (data not shown), all individuals sampled 
from this region were analyzed to estimate the genetic charac- 
teristics of Kangaroo Island's Victorian founder population. 
Genetic analysis also revealed similarity between individu- 
als of Tasmania's northwestern catchments (data not shown) 
and, consequently, all individuals sampled from this region 
were used to estimate the genetic characteristics of Kangaroo 
Island's Tasmanian founder population. 

For the microsatellite data, fstat version 2.9.3 (Goudet 
2001) was used to calculate allelic richness averaged over 
loci and Weir and Cockerham's (1984) measure of F\ S . Sig- 
nificant differences in allelic richness were calculated by a 
Wilcoxon signed-rank test (Wilcoxon 1945). Observed (H 0 ) 
and expected heterozygosities (H E ) were estimated and devi- 
ations from Hardy-Weinberg equilibrium (HWE) were de- 
termined by exact tests and permutation in arlequin version 
3.11 (Excoffier et al. 2005). A global estimate and population 
pairwise estimates of F ST were calculated in fstat (1500 per- 
mutations). A factorial correspondence was implemented in 
genetix version 4.05 (Belkhir et al. 2004) to visualize dif- 
ferences in genetic diversity between populations. The two 
factors that explained the majority of the variation in multi- 
locus genotypes were plotted. 

Structure version 2.3.3 (Pritchard et al. 2000) was run to 
determine the relationship between individuals. Using only 
genetic data, Structure identifies the number of distinct 
clusters, assigns individuals to clusters, and identifies ad- 
mixed individuals. To determine the number of populations 
(K) within the data set, five independent simulations of each 

(1—5) were run with 100,000 burn-in iterations and 500,000 
data iterations. The most likely value of K is the one that max- 
imizes the log-likelihood of obtaining the observed sample 
of multilocus genotypes (Pritchard et al. 2000). We used the 
method of Evanno et al. (2005) to estimate the true K. Struc- 
turama (Huelsenbeck et al. 2011) was also used to infer the 
number of populations. This program differs slightly from 
structure in that it allows you to place a prior distribution 
on K and estimate the number of populations directly. 

To determine the contribution of founders on Kangaroo 
Island, a genetic assignment test was performed using the pro- 
gram GeneClass2 (Piry et al. 2004). Tasmanian and Victorian 
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individuals were partitioned to create two independent refer- 
ence populations. Kangaroo Island individuals were assigned 
to a reference population based on the Bayesian methods of 
Rannala and Mountain (1997). The probability of correctly 
assigning individuals to their population of origin was calcu- 
lated from 1,000,000 Monte Carlo simulations with a Type I 
error of 0.01. High values indicate a good fit to the assign- 
ment population, while values less than about 0.05 suggest a 
poor fit to the population. 

Sequenced mitochondrial products were aligned using 
mega version 4.0 (Tamura et al. 2007) and edited manually. 
Levels of genetic differentiation between the haplotypes were 
calculated using pairwise genetic distances under a Kimura- 
2-parameter model in mega version 4.0 (Tamura et al. 2007). 
A haplotype network was constructed in tcs 1.21 (Clement 
et al. 2000) to visualize base pair changes and relationships 
between haplotypes. 

Population size modeling — Kangaroo island 

Intensive platypus mark-release-recapture studies were car- 
ried out along Rocky River in Kangaroo Island. The surveyed 
area consists of relatively high quality habitat (R. Ellis, un- 
publ. data) and is likely to contain proportionally more platy- 
puses than regions upstream or downstream. Consequently, 
despite the surveyed area constituting approximately one- 
third of the platypus' known distribution on the island, it 
is estimated to contain approximately half the population 
(R. Ellis, unpubl. data). Surveys were carried out roughly 
every month from 1 March 1998 until 19 November 1999 with 
a total of 92 capture events recorded across 44 unique adult 
individuals. For analytical purposes, each sampling month 
was allocated to a capture occasion. Sex, body weight, head- 
to-tail length, bill width, bill length, and shield measurements 
were also recorded for each captured individual. These mea- 
surements were normalized prior to model fitting. Data for 
males and females were analyzed separately. For the female 
data there were no records between March 1998 to July 1998, 
May 1999 to July 1999, and for October 1999. To avoid data 
sparseness, observations from August 1998 and September 
1998 were pooled. This gave 40 observations of 20 unique 
female platypuses across 1 1 capture occasions. For the male 
data, there were no records for May 1998, November 1998, 
February 1999, June 1999, and October 1999. We pooled 
observations for male data for the months March/ April for 
both years 1998 and 1999. This gave 52 observations of 
24 unique male platypuses across 14 capture occasions. Indi- 
vidual capture histories were then constructed and covariates 
were taken as the first recorded observation and assumed 
constant across each capture occasion. 

The population size for platypus within the surveyed area 
was denoted as N. An estimate of the population size was 
denoted by N and their standard errors by SE ( N), which can 



be used to obtain 95% confidence intervals (CIs). Each data 
set was tested for closure because different models are used 
when the population is assumed closed or open. For closed 
populations, the population is assumed to remain constant 
throughout sampling. For open populations, the population 
is susceptible to births, deaths, immigrations, and emigration 
and if individuals leave the population, they do so perma- 
nently. To test for closure we consider the tests of Stanley and 
Richards (1999) where we define the Stanley and Burnham 
(1999) test as Test 1 and the Otis et al. (1978) test as Test 2. In 
these tests, the null hypothesis assumes closure, so significant 
(a < 0.05) P values suggest the population is not closed. 

In closed population models, only the population size and 
capture probabilities are estimated. Conditional likelihood 
models were used (Huggins 1989, 1991) where capture prob- 
abilities can permit heterogeneity between individuals (h), 
behavioral response (b), or time (t) (Otis et al. 1978). That 
is, capture probabilities are modeled linearly as a function 
of: individual covariates, such as body weight or gender for 
heterogeneity models, time-effects for time-dependent mod- 
els, and trap behavior effects for behavioral response models. 
These models are denoted by M(„), for example, model M(,b) 
allows for capture probabilities to depend on time and behav- 
ioral response. Also, model M (0 ) assumes all individuals have 
the same probability of capture. There are eight combinations 
of these models in total. 

For open population models, all individuals are assumed to 
have the same probability of capture. In addition to estimat- 
ing the population size, survival and recapture probabilities 
and new arrivals to the population are also estimated. The 
classical Jolly-Seber (JS) model (Seber 1982; Pollock et al. 
1990) can be used, however, we use the more modern ap- 
proach of Schwarz and Arnason (1996), which generalizes 
the way birth is modeled in the usual JS model. This allows 
an open population estimate to be obtained. 

For both models, the estimates are obtained by using cap- 
ture histories and maximum likelihood estimation (section 
1.3 in Amstrup et al. [2005]). The Akaike information crite- 
rion (AIC) statistic (Burnham and Anderson 1999) was used 
for model selection. 

The average effective population size over time was esti- 
mated from genetic data following the equation of Culver 
et al. (2008). 



1 



1 

2N e 



where H t = expected heterozygosity in generation t , H 0 = 
initial heterozygosity, and N e = the effective population size. 

The effective population size on Kangaroo Island was cal- 
culated according to the equation N e /N = 0.1 (Frankham 
1 995b ). A simple simulation was undertaken for Kangaroo Is- 
land's population to predicted changes in heterozygosity and 
allelic richness over time. The simulation was undertaken in 
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Table 1. Thirteen microsatellite loci screened, primer sequences, annealing temperatures, and the number of alleles found in the contemporary 
populations of upper Yarra Victoria, northwestern Tasmania; Kangaroo Island, South Australia; and King Island, Tasmania. Unique alleles are shown in 
brackets with alleles unique to northwestern Tasmania in italics and alleles unique to Victoria's upper Yarra population in bold text. 



Locus 



Repeat motif 



Primer sequence 5'-3' 



Upper Yarra Northwestern Kangaroo King 
7~ a Victoria Tasmania Island Island 

(n = 60) (n = 19) (n = 12) (n = 21) 



OA 1 .3 (GA) 5 N (AG) 5 N 20 (AG),, 

OA 3.2 (GA) S (AT) S N,„ (TG), 2 

OA 4.5 (AC)„ 

OA 5.1 (GT) 14 

OA 6.2 (TC) 22 

OA 7.3 (TG) 12 

OA 10.5 (TC) 9 N (CA),„ 

OA 11.9 (GA) 5 (GT) 10 

OA 12.6 (CA) 14 

OA 14.3 (CA) 14 (CACAC) 2 

(AC) 5 (CA) 13 

OA 17.6 (GA) 10 

OA 18.5 (TG),, N 3 (GC) 6 

OA 20. 12 (TC), 5 (CT) 17 



F: GACCTCTTTGCCACTTTGCTA 59.4 10(6) 5(7) 

R: GGATTAGAACCCACGATCTGTT 

F: GCCCTATGTACCTTGAATATAA 48.7 6 (3) 3 

R; ACAGTTGGTGGACTTGATTC 

F: ACGCCCCACCCGTTCCCTTTC 61.4 6(2) 7(3) 

R: ATCCATTCGCCGATCTCCTGTGC 

F: CTTGGAAAGCATACACAGATG 52.0 4 (2) 2 

R; GAAATTGTTGGACTATGGGTAT 

F: TAGGGTGGTTTGAAAGGTTTTG 56.1 18(14) 9(5) 

R; AGACAGCCGTAGGAGCACTAAA 

F: AATCTGAAAAGGCAACAATCT 48.4 12(8) 5(7) 

R; GGGCTTATCATTTGTCCTCTA 

F: GCTCTGATGGCTAATACTGCTA 50.3 3 3 

R; ATCCCTTCCCTCTCCATTATTA 

F: GGTCAAAGAGTCCCAGAATGAC 60.1 1(1) 2(2) 
R: GAGACAGGAAACTTGGCATAGG 

F: GATCTCCCACTACCGACAGTTT 56.1 9(7) 3(7) 

R; CAGGGTGGAATGATTACAGAAA 

F: GAAGGAGGAGGAGAGGTTGACA 54.0 14(10) 5(7) 
R: TTCAGCGACTTTTCTGTTCCATAG 

F: GTAACTTCTCACGGGGCAACTT 53.5 8 (6) 3 (7) 

R: G G C ATTTTATTTTC TCGCCTCTA 

F: TTGTCTATATTCTTGGAAGGGCTC 60.0 10(6) 5(7) 
R; ATTGCAGGTAAAGTGAAGGGAA 

F: GTTCCCTTGAGGACGGAGA 60.1 16(7) 11(2) 

R: CAGTGGGCTTTTCCATTCATA 



2(1) 

2 

2 

3(1)(1) 
5(5) 
3(1) 
1 

1 (D 
5(4) 
3(2) 
2 

3(1) 
5(3) 



1 
2 
1 

1 (D 
1 (D 
2(1) 
1 

1 (D 
3(1) 
1 

2 

1 (D 
2 



7"a, annealing temperature; number of alleles observed (number of unique Tasmanian alleles are given in parenthesis). 

* Allele not present in either upper Yarra or northwestern Tasmanian individuals sampled — broader sampling has detected this allele in other contem- 
porary mainland populations including <50 km from the upper Yarra population (data not shown). 



PopTools (Hood 2002) using random sampling of allele fre- 
quency data each generation for a given N e . The simulation 
assumes random mating, nonoveiiapping generations, and a 
constant effective population size across time. 

Results 

Genetic diversity 

A total of 1 36 alleles were detected across 1 3 polymorphic mi- 
crosatellite loci. All loci were polymorphic in Tasmania while 
Victoria had one monomorphic locus (OA 11.9), Kangaroo 
Island had two (OA 11.9 and OA 10.5), and King Island had 
eight monomorphic loci (Table 1). The majority of all alleles 
identified (>85%) were present in the Victorian population. 
There was a clear genetic distinction between individuals of 
Victoria and Tasmania with unique alleles observed at 12 
loci. Victoria's population contained 72 unique alleles and 
Tasmania's population contained 18 unique alleles (Table 1). 



No unique alleles were present in the King Island samples; 
all alleles were present in Victoria's upper Yarra population 
and incorporated several alleles unique to Victoria. Simi- 
larly, nearly all the alleles identified from Kangaroo Island 
were present in Victoria's upper Yarra individuals and in- 
corporated several alleles unique to Victoria. One allele was 
detected on Kangaroo Island unique to this population. Nei- 
ther Kangaroo Island nor King Island contained any of the 
alleles identified as being unique to Tasmania. 

Genetic diversity estimates were similar in the Tasmanian 
and Victorian populations; expected heterozygosities were 
0.606 in Tasmania and 0.597 in Victoria. Allelic richness es- 
timates differed significantly between Victoria and Tasmania 
(P < 0.01) with values of 4. 34 and 5.79, respectively (Table 2). 
Genetic diversity estimates of the contemporary Kangaroo Is- 
land population (i.e., sampled 2008-2009) revealed a lower 
level of allelic richness (r = 2.85) and expected heterozygos- 
ity (H E = 0.419) compared with Tasmanian and Victorian 
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Table 2. Ornithorhynchus anatinus population genetics statistics for individuals sampled from the upper Yarra Victoria, northwestern Tasmania, 
King Island Tasmania, and South Australia. Samples from South Australia have been divided into categories to calculate genetic diversity in each of 
two groups: the most recent Kangaroo Island population (sampled from 2008 to 2009) and all Kangaroo Island and Warrawong Sanctuary individuals 
(sampled from 1991 to 2009). All individuals were genotyped at 13 microsatellite loci. 



Region 


n 


a 


r 


Ho 


H E 




HWE P-value 


rA 






Upper Yarra, Victoria 


60 


9.00 


5.76 


0.549 


0.597 


0.08 


0.646 


5.077 


(±3.499) 


Northwestern Tasmania 


19 


4.85 


4.34 


0.550 


0.606 


0.095 


0.606 


1.154 


(i 


1 .463) 


King Island, Tasmania 


21 


1.46 


1.32 


0.026 


0.032 


0.205 


0.083 


0.462 


(±0.660) 


Kangaroo Island, South Australia, 2008-2009 


12 


2.85 


2.85 


0.423 


0.419 


-0.01 


0.495 


0.692 


(±0.947) 


Kangaroo IslandAA/arrawong Sanctuary, South Australia 


25 


3.54 


3.02 


0.395 


0.431 


0.086 


0.467 


1.077 


(± 


1.256) 



n, the number of individuals genotyped for each population; a, mean number of alleles; r, allelic richness; H 0 , observed heterozygosity; H E , expected 
heterozygosity; FIS, multilocus estimates; HWE, Hardy-Weinberg equilibrium P-values. 

*Significance after corrections for multiple comparisons. Locus; r/\, mean number of rare alleles (frequency < 0.05) per locus. 



populations. King Island exhibited extremely low genetic di- 
versity, with an allelic richness of 1.32 and expected het- 
erozygosity of 0.032 (Table 2). Inbreeding levels were also 
high within the King Island population (F K = 0.205, P < 
0.01) but the population was not out of HWE, which is likely 
due to low sample size and extremely low numbers of al- 
leles. In comparison to Victoria, King Island and Kangaroo 
Island contained few rare alleles (mean 0.462 and 0.692, re- 
spectively). Numbers of rare alleles in Tasmania were also 
comparatively low (mean 1.154). 

Significant genetic differentiation was present between the 
Tasmanian and Victorian populations with an P ST of 0.254 
(95% CI, 0.144-0.394). King Island was significantly differ- 
entiated from these two populations with an F ST of 0.396 
with Victoria (95% CI, 0.297-0.478) and 0.635 with Tas- 
mania (95% CI, 0.554-0.715). South Australian platypuses 
(Kangaroo Island and Warrawong individuals) were only 
slightly differentiated from Victorian samples (F S t = 0.079, 
95% CI, 0.043-0.114) but highly differentiated from Tas- 
manian samples (P S t = 0.366, 95% CI, 0.266-0.476). The 
factorial correspondence analysis provides a visualization of 



these differences (Fig. 2), showing the similarity between 
individuals of Victoria and Kangaroo Island and the clear 
distinction between individuals of Tasmania, Victoria, and 
King Island. 

Structure analysis indicated the data set was comprised 
of three populations (K = 3) and this was supported by 
Structurama. King Island individuals were allocated to a 
distinct population group, as were Tasmanian individuals. 
The third population comprised all Kangaroo Island, War- 
rawong, and Victorian individuals (Fig. 3). Assignment tests 
assigned all Kangaroo Island individuals to Victoria with a 
mean probability of 0.86 (range: 0.46-1.00). Tasmania was 
excluded as the population of origin for all Kangaroo Island 
individuals. 

mtDNA sequencing produced a clear separation between 
Victorian and Tasmanian individuals with 1 1 fixed nucleotide 
differences at COII and 16 fixed differences at cytb. Only one 
COII and one cytb haplotype was present in individuals from 
Victoria (Fig. 4). The same haplotype was observed in all 
Kangaroo Island individuals analyzed. Tasmania contained 
more mtDNA variation with four COII haplotypes and six 



Figure 2. Factorial correspondence analysis 
depicting the microsatellite genetic difference 
between Ornithorhynchus anatinus individuals 
of four locations: the upper Yarra, Victoria; 
northwestern Tasmania; King Island, Tasmania; 
and Kangaroo Island and Warrawong 
Sanctuary, South Australia. Differences are 
represented across two factors with factor 1 
representing 51 .07% of the variation and 
factor 2 representing 37.79% of the variation. 
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Figure 3. Summary plot of the estimated membership coefficient for each Ornithorhynchus anatinus individual across three population clusters (K 
= 3). Each individual is represented by a single vertical line broken into three shaded segments. Segment represents the proportional membership to 
each of the three population clusters. Individuals are grouped into the regions from which they were sampled: upper Yarra, Victoria; northwestern 
Tasmania; Kangaroo Island and Warrawong Sanctuary, South Australia; and King Island, Tasmania. 



Figure 4. Mitochondrial DNA network 
analysis showing the relationship between 
haplotypes for mtDNA gene (a) COII and (b) 
cytb. Dots on branches indicate nucleotide 
changes. Gray haplotypes indicate presence in 
Victoria's upper Yarra while white haplotypes 
indicate presence in northwestern Tasmania. 
Black stripes identify haplotypes present in 
Kangaroo Island and Warrawong individuals 
and black dots indicate haplotypes present in 
King Island, Tasmania. The number of 
individuals revealing each COII haplotype are 
as follows: A = SO; B = 32; C = 3; D = 1 ; E = 
1 and for cytb: A=11;B=11;C = 4;D = 3; 
E = 1; F = 1; and G = 1. 




cytb haplotypes (Fig. 4). King Island individuals' contained 
only one haplotype at each mitochondrial gene, both of which 
were present in Tasmania. Kimura-2-parameter distances be- 
tween COII Victorian haplotype A and Tasmanian haplotypes 
ranged between 0.0257 and 0.0277, while distances between 
cytb Victorian haplotype A and Tasmanian haplotypes ranged 
from 0.0236 and 0.0293. Genetic differentiation within Tas- 
manian samples was low with COII haplotype (B, C, D, E) 
differences ranging between 0.0019 and 0.0098, while cytb 
haplotype (B, C, D, E, F, G) differences ranged from 0.0014 
to 0.0124. 



Population size estimates — Kangaroo island 

Each data set (females only and males only) was tested 
for closure (Table 3). The female data did not indicate 



Table 3. Closure tests for the female and male data. Low P values 
indicate lack of closure. 



Data 


Test 1 , P-value 


Test 2, P-value 


Conclusion 


Male 


0.02 


0.49 


Possibly open 


Female 


0.08 


0.40 


Closed 
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Table 4. Various parametric models fitted to female and male Kangaroo Island platypus data with the body weight covariate for M (h) -type models. 
Note that M (h 2, represents a quadratic model. 



Data 


Model 


M m 


M m 


M m 


M„ 






M w 






Female 


AIC 


221.18 


204.40 


219.55 


220.52 


203.41 


204.75 


202.59 


218.56 


201.46 




N 


20.75 


21.73 


20.90 


24.21 


24.29 


24.74 


21.83 


24.18 


24.27 




SE(W) 


1.42 


2.17 


1.68 


2.92 


2.97 


3.64 


2.28 


2.90 


2.94 


Male 


AIC 


299.38 


286.09 


299.55 


297.71 


285.88 


284.74 


286.22 


297.76 


285.87 




W 


25.98 


25.58 


25.75 


28.15 


28.43 


32.22 


25.27 


27.42 


27.68 




SE(/V) 


3.05 


1.94 


2.93 


2.95 


3.07 


6.25 


1.63 


2.41 


2.52 



lack of closure, thus, we assume the female population is 
closed. For the male data, the two tests reported differ- 
ent conclusions; therefore, open and closed models were 
considered. 

For the female data, all eight closed population models 
of Huggins (1991) were fitted using the body weight covari- 
ate for M(h) type models. Other covariates were not used to 
avoid collinearity and body weight measurements were nor- 
malized prior to model fitting. Based on the AIC, model M (0 ) 
provided the best fit but models M(b) and M(h) may also be 
appropriate (Table 4). In order to compare the models fit- 
ted to female and male data (see below), a further analysis 
of model M(h) was carried out by considering a quadratic 
relationship between the covariate and capture probability, 
however, this did not improve the model fit. Fitted capture 
probabilities for Model M(h)were plotted against body weight 
with CIs. The relationship was very flat and close to a constant 
model, which may explain why model M( 0 ) provided a better 
fit. According to model M( 0 ), the population size of females 
within the sampled area on Kangaroo Island was estimated 
as N(Female) = 24.27 (18.39, 30.15) (Table 4). 

For the male data, both closed and open population 
models were fitted. Initially, the open population model 
described above using program MARK (White and Burn- 
ham 1999) was fitted. The open population estimate was 
JV(Male) = 30.86 (26.18, 45.61), however, identifiability is- 
sues were encountered for most model parameters due to 
sparseness of data. The results for closed population models 
that also use a normalized body weight covariate are shown 
in Table 4. The best selected models were M( 0 ) and M(h) . A 
quadratic model was fitted, which gave a lower AIC and indi- 
cated some nonlinearity. If we consider a closed population, 
the estimated number of male platypuses within the surveyed 
area was N(Male) = 32.22 (19.72, 44.72). As male platy- 
puses are typically more mobile than females (Grant 2007), 
it may be more appropriate to consider the male population 
open. In any case, the estimated population sizes for closed 
and open population models are similar at approximately 
31-32 individuals. 

In total, a population of approximately 55 individuals was 
estimated within the study area. As this study area is assumed 



to contain approximately half the platypus population of 
the island, the total population of the island is expected to 
be roughly double this figure (i.e., ~1 10 individuals). This 
equates to an effective population size of around 1 1 individ- 
uals {NJN =0.1). 

To estimate the effective population size from genetic data, 
we use the formulae of Culver et al. (2008). The genetic diver- 
sity present in today's Kangaroo Island population is 0.419 
(H t = 0.419). We expect the initial level of heterozygosity in 
the Kangaroo Island population to be similar to that present 
in the founding population in the upper Yarra River today 
(H 0 = 0.597). Given platypuses are capable of breeding in 
their second year and have been recorded surviving up to 
21 years of age (Grant 2004), the generation time is assumed 
to be ~10 years (although this is likely to contain some de- 
gree of error). Approximately seven generations (f = 7) have 
passed since Victorian platypuses were first introduced to 
Kangaroo Island. This equates to an effective population size 
of 10.14 or ~ 10 individuals, a value similar to that produced 
by population size models. 

Using a simple simulation with effective population size 
(N e ) of 1 1 , genetic diversity is predicted to continue to decline 
within Kangaroo Island's platypus population through time 
(Fig. 5). Within 50 generations, expected heterozygosity and 
average allele numbers are anticipated to be similar to those 
observed in King Island's platypus population. 

Discussion 

Levels of genetic diversity in platypus populations from Tas- 
mania and Victoria are similar to those observed in mainland 
populations of other Australian mammalian species (Taylor 
et al. 1994; Pope et al. 2000; Bowyer et al. 2002; Eldridge et al. 
2004). The sampled region of Tasmania does, however, pos- 
sess a reduction in allelic richness and lower frequency of rare 
alleles in comparison to Victoria's upper Yarra population 
but mitochondrial DNA diversity remains high. On Kanga- 
roo Island, levels of genetic diversity are reduced with ~30% 
decrease in expected heterozygosity and ~50% decrease in 
allelic richness in comparison to Victoria. The genetic diver- 
sity of King Island is extremely low with a high proportion of 
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Figure 5. Simulated loss of genetic diversity 
over time (in generations) of Ornithorhynchus 
anatinus individuals on Kangaroo Island, 
South Australia (H E , average heterozygosity 
and A, average allele number). The population 
was assumed to have commenced with levels 
of genetic diversity similar to that found in the 
present day upper Yarra population in Victoria. 
One thousand Monte Carlo simulations were 
run for 13 microsatellite loci assuming random 
mating, nonoverlapping generations and A/ e = 
1 1 . Solid lines indicate means and dotted lines 
represent 95% confidence intervals. The 
current genetic diversity of the Kangaroo 
Island population is indicated by a black 
diamond. 




monomorphic loci, ~70% reduction in allelic richness, and 
> 90% reduction in expected heterozygosity in comparison to 
Victoria. These microsatellite genetic diversity estimates for 
King Island platypuses are among the lowest ever recorded 
for a naturally outbreeding vertebrate population. 

Island populations frequently have reduced genetic diver- 
sity in comparison to mainland populations (Frankham 1997, 
1998). The rate of loss of genetic diversity is dependent on 
the effective population size (N e ) and the number of genera- 
tions over which the population has been isolated (Frankham 
1997). Consequently, reductions in genetic diversity are of- 
ten more pronounced in endemic rather than relatively re- 
cently introduced or translocated populations (e.g., Cardoso 
et al. [2009]). King Island has been isolated from mainland 
Australia for close to 12,000 years and isolated from Tasma- 
nian for ~ 10,000 years (Blom 1988). This long-term isolation 
combined with small population size has likely contributed 
to the severe reduction in genetic diversity found in this 
O. anatinus population. While no gross morphological 
changes have been detected (i.e., length, weight, external 
appearance), additional phenotypic traits have not been as- 



sessed within this population and therefore it is not known if 
the severely low levels of genetic diversity are impacting pop- 
ulation fitness (inbreeding depression). In any case, the low 
levels of genetic diversity suggest that this population is ex- 
posed to an increased risk of extinction (Markert et al. 2010), 
and if it is challenged by environmental stress and/or parasites 
and disease, platypuses may become extinct on King Island. 

Determining the relationship between King Island's platy- 
puses with those of Tasmania or Victoria is difficult because 
contact with the major land masses to the north and south 
was lost ~10,000 years ago. Although King Island retains a 
mtDNA haplotype present in Tasmania, it also contains sev- 
eral nuclear alleles that are unique to today's Victorian popu- 
lation. Thousands of years of independent evolution will have 
seen allele frequencies and distributions alter significantly on 
both island and mainland populations. For the small King 
Island population, this extended period of isolation has also 
lead to extremely low genetic diversity, low allele number, 
and consequently, its identification as a distinct population. 

The greater proportion of heterozygosity retained on Kan- 
garoo Island (~70%) is likely to reflect this population's 
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relatively recent isolation, around 70 years ago. The present- 
day population on Kangaroo Island, however, is likely to 
be limited in its population size. From a small number of 
founders, the Kangaroo Island population has undergone 
population growth, attaining an estimated size of ~ 1 10 indi- 
viduals by 1998/1999. An extremely low proportion of juve- 
niles encountered during trapping periods ( <5% in compar- 
ison to nearly 33% on the mainland [Grant 2004] ) appears 
to indicate that the population may have reached carrying 
capacity within its current distribution along the Rocky and 
Breakneck Rivers. Despite records of platypuses traversing 
overland on the island (Ellis 2000), individuals have yet to 
be encountered in waterways outside of the two they were 
initially introduced into. There maybe potential for the pop- 
ulation to expand their range into other waterways on the is- 
land and thus increase the total population size. Within their 
current distribution, however, the population is unlikely to 
grow beyond 110 individuals. The small number of founders, 
genetic isolation, and sustained small population size will 
continue to impact on the genetic diversity of this island 
population. Both microsatellite and mitochondrial data only 
provide support for the contribution of Victorian founders 
suggesting the initial Tasmanian founders were unsuccessful. 
This limits the founder population to a maximum of 16 Vic- 
torian individuals and potentially fewer effective founders 
(those that have contributed genetically to subsequent 
generations). 

Today's Kangaroo Island platypus population remains 
small (~JV = 110) and the number of individuals contribut- 
ing to subsequent generations is likely to be even smaller. 
Estimates of the effective population size of Kangaroo Island 
obtained from both population size modeling and genetic 
data were similar, indicating an N e of approximately 10— 
1 1 individuals. This is extremely low and potentially unsuit- 
able for long-term population persistence. An N e of at least 
500 (Franklin 1980), but more likely greater than 1000 (Willi 
et al. 2006; but see [Lande 1995]) has been proposed as a 
minimum to maintain the evolutionary potential of a pop- 
ulation. At N e < 50, inbreeding is likely to be high and can 
limit a population's ability to adapt even in the short-term 
(Franklin 1980). Recently bottlenecked populations are likely 
to experience a reduction in allele number, particularly due 
to the loss of rare alleles, but can still maintain reasonable lev- 
els of heterozygosity in the short-term (Luikart et al. 1998). 
This is in line with the findings on Kangaroo Island where 
allelic richness (r = 2.85) and the proportion of rare alleles 
are low (rA = 0.692) compared to Victoria, whereas heterozy- 
gosity levels are reduced but not to the same extent (H E = 
0.419). Isolation since 1946 in combination with a small ef- 
fective population size is likely to have elevated the impact of 
genetic drift, contributing to the observed decrease in allele 
number (particularly the loss of rare alleles) and decrease in 
heterozygosity. 



With a small population size and continued isolation of 
the Kangaroo Island population, genetic diversity is likely 
to continue to diminish through time (Fig. 5). Without in- 
tervention, it is likely that genetic diversity will continue to 
decline and be in a similar state to that observed on King 
Island within 50 generations (~200-500 years according to 
generation times based on the average age of breeding females 
and longevity estimates in Grant [2004]). 

Long-term isolation and low N e are also likely to have led 
to inbreeding on these island populations. Islands typically 
exhibit increased levels of inbreeding with the effective in- 
breeding (F e ) estimated from equation 1 from (Frankham 
1998), 



where H K is the heterozygosity of the island population and 
H M is the heterozygosity of the mainland population. For 
O. anatinus on Kangaroo Island, the effective inbreeding co- 
efficient is large (_F e = 0.30) and for the King Island popula- 
tion, inbreeding is severe (F e = 0.95). 

Low levels of genetic diversity and high levels of inbreeding 
are likely to have important consequences for the long-term 
survival of island populations. Both measures correlate with 
population persistence as well as various fitness indicators 
including reproductive success, survival, and parasite resis- 
tance (Ralls et al. 1988; Westemeier et al. 1998; Coltman et al. 
1999; Crnokrak and Roff 1999; Seymour et al. 2001; Brook 
et al. 2002; Reed and Frankham 2003). Despite the occurrence 
of long-term inbreeding on King Island, the genetic load is 
still likely to be high (Frankham 1995a; Ballou 1997; Eldridge 
et al. 1999). Indeed, King Island individuals harbor low di- 
versity in an important immune response gene family, the 
major histocompatability complex (MHC) (Lillie 2009). The 
occurrence of the fungal mucormycosis on the Tasmanian 
mainland (Gust and Griffiths 2009) is of particular concern 
and could have potentially devastating consequences should 
it emerge within the island population. High levels of ge- 
netic diversity are required to maintain adaptive potential 
and minimise the risk of extinction (Reed and Frankham 
2003). 

Conclusion 

Islands can act as an important reservoir for Australian 
wildlife. They are often free from introduced predators or 
competitors. There are numerous examples of species persist- 
ing on islands despite extinction on the mainland (Van Dyck 
and Strahan 2008). Island populations can provide a valuable 
source of individuals for conservation programs, including 
captive breeding, reintroductions, or translocations. How- 
ever, island populations also present many problems that 
must be effectively managed if they are to be employed in 
conservation efforts. 
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Currently, genetic diversity in King Island O. anatinus is 
severely depauperate and the population is likely to be suffer- 
ing from reduced fitness, reduced evolutionary potential, and 
an increased risk of extinction. Consequently, their long-term 
viability is likely to be under threat. Levels of genetic diversity 
on Kangaroo Island are reduced in comparison to mainland 
populations and are predicted to decrease to levels as low as 
King Island within approximately 50 generations (Fig. 5). To 
maintain adaptive potential and minimise the risk of extinc- 
tion (Reed and Frankham 2003), levels of genetic diversity 
need to be maintained (in the case of Kangaroo Island) or 
ideally, increased. If additional suitable platypus habitat can 
be found on Kangaroo Island, increasing the total popula- 
tion size through population range expansion can slow the 
loss of genetic diversity. Human-mediated migration pro- 
vides an alternative means to mitigate the negative effects 
of lowered genetic diversity with the transfer of individuals 
from suitable source populations. The contribution of just 
one effective migrant per generation is presumed sufficient 
to alleviate the effects of drift and reduce inbreeding depres- 
sion while maintaining local adaptation (Wang 2004). The 
long-term persistence of these populations will ultimately 
depend on adequate levels of genetic diversity (Bouzat 2010) 
and the addition of new genetic material provides a means to 
achieve that goal. 
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